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Q_i We consider a model of scalar field with non minimal kinetic and Gauss Bonnet 

couplings as a source of dark energy. Based on asymptotic limits of the general- 
ized Friedmann equation, we impose restrictions on the kinetic an Gauss-Bonnet 
couplings. This restrictions considerable simplify the equations, allowing for ex- 
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(f-) act solutions unifying early time matter dominance with transitions to late time 

quintessence and phantom phases. The stability of the solutions in absence of 
Q^ matter has been studied. 
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1 Introduction 

One of the most important challenges for the current theory of cosmology is the expla- 
nation of the late-time accelerated expansion of the universe. The cosmic acceleration 
has been supported by different observational data such as the supernovae type la 
[I], 0, [2], [1], cosmic microwave background anisotropy and large scale galaxy 
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surveys [6]. The interpretation of astrophysical observations indicates that this accel- 
erated expansion is due to some kind of negative-pressure form of matter known as 
dark energy (DE). The combined analysis of cosmological observations also suggests 
that the universe is spatially flat, and consists of about ~ 1/3 of dark matter, and 
~ 2/3 of homogeneously distributed dark energy with negative pressure. The simplest 
candidate for dark energy is the cosmological constant, conventionally associated with 
the energy of the vacuum. However the cosmological constant presents the known 
severe energy scale problem [7], [8]. Among the different alternatives to explain the 
dark energy, the models involving a scalar field play an important role as they provide 
a dynamical behavior of the equation of state, which is favored by astrophysical data 
(for a review see [9j Ell EJ EE] ). A widely explored are scalar-tensor models, which 
contain a direct coupling of the scalar field to the curvature, providing in principle a 
mechanism to evade the coincidence problem and naturally allowing (in some cases) 
the crossing of the phantom barrier [T3], [2], [T5]. From a pure geometrical point 
of view, the modified gravity theories, which are generalizations of the general rela- 
tivity, have been widely considered to describe the early-time inflation and late-time 
acceleration, without the introduction of any other dark component, and represent 
an important alternative to explain the dark energy (for review see |16j). 
In the present work we consider a string and higher-dimensional gravity inspired scalar 
field model, with two types of couplings: kinetic coupling to curvatures and Gauss 
Bonnet (GB) coupling to the scalar field. The GB term is topologically invariant in 
four dimensions and does not contribute to the equations of motion. Nevertheless 
it affects cosmological dynamics when it is coupled to a dynamically evolving scalar 
field through arbitrary function of the field, giving rise to second order differential 
equations of motion, preserving the theory ghost free. Therefore, the coupled GB 
term seems as a natural generalization of the scalar field with non-minimal kinetic 
coupling to curvature [T71 HE] • 

Some late time cosmological aspects of scalar field model with derivative couplings to 
curvature have been considered in [TH], [13 [EE1 [20] , [21], [22], [23]- On the other hand, 
the GB invariant coupled to scalar field have been extensively studied. In p3], [25] 
the GB correction was considered to study the dynamics of dark energy, where it was 
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found that quintessence or phantom phase may occur in the late time universe. The 
evolution of perturbations in GB dark energy at large and small scales has been cal- 
culated in [25], [26J. Accelerating cosmologies with GB correction in four and higher 
dimensions have been discussed in [27], [28], [29] . The modified GB theory applied 
to dark energy have been suggested in [30] , and different aspects of the modified GB 
model applied to late time acceleration, have been considered among others, in [3Tj . 
[32], [33], [31], [35]. 

All these studies demonstrate that it is quite plausible that the scalar-tensor cou- 
plings predicted by a fundamental theory like string theory, may become important 
at current low-curvature universe (see [36J and [H] for review). 
The most general second-order ghost-free scalar-tensor Lagrangian with couplings to 
curvature, can be originated from toroidal compactification of 4 + A dimensional La- 
grangian of pure gravity, as shown in [37]. In this compactification the scalar field 
plays the role of the overall size of the A-torus, and the couplings appear as ex- 
ponentials of the scalar field. This general Lagrangian also appears in the next to 
leading order corrections in the a' expansion of the string theory [38], [39], [ID]. In 
the present study we consider those terms in this ghost-free scalar-tensor Lagrangian 
that are coupled to the curvature, namely kinetic terms coupled to curvature and 
scalar field coupled to GB invariant, where we use more general couplings in order 
to find appropriate cosmological scenarios that satisfy the current astrophysical ob- 
servations. Some late-time cosmological solutions have been studied in [JT], and Big 
Rip and Little Rip solutions have been reconstructed in [42]. The main purpose of 
this work is to find exact solutions that are consistent with the large amount of exist- 
ing data supporting the late-time accelerated expansion of the universe. Even in the 
Friedmann-Robertson- Walker background the cosmological equations are very diffi- 
cult to integrate in an exact way, in part due to their non linearity and the variety of 
conditions that should be satisfied in order to cover the different evolutionary stages. 
Nevertheless, introducing restrictions on the interaction terms that are appropriate 
and consistent with early and late time cosmology, we can find suitable exact solutions 
describing early and late time asymptotic behavior. It was shown that the model with 
string (exponential) couplings in absence of matter, leads to power-law solution. In 
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section II we introduce the model and give the general equations, which are then ex- 
panded on the FRW metric. In section III we consider two type of restrictions on the 
GB and Kinetic couplings in the case of scalar field dominance and found the respec- 
tive solutions describing power-law, and more general viable late-time cosmologies. 
The issue of stability of these solutions is also addressed. In section IV we consider 
the matter contribution and study exact solutions unifying early time decelerated 
with late time accelerated expansion, including phantom phase. Concluding remarks 
are given in section V. 



2 Field Equations 

Let us start with the following action containing scalar field with kinetic couplings to 
curvature and the Gauss Bonnet coupling 

S = [ d 4 x^-g\-^—R - l-d^dy + F l {<j>)G tlv d^d v ( j> - V^) - F 2 (0)d + S m 
J Llb7rG 2 J 

(2.1) 

where G^ = — ^g^R, Q is the 4-dimensional GB invariant Q = R 2 — AR^R^ + 
R^vpcrR^^ and S m is the matter action. The coupling Fi(<p) has dimension of 
(length) 2 , and the coupling F 2 (<f)) is dimensionless. The GB coupling does not make 
contributions higher than second order (in the metric) to the equations of motion, 
and therefore does not introduce ghost terms into the theory. Hence, the equations 
derived from this action contain only second derivatives of the metric and the scalar 
field, avoiding problems with higher order derivatives [33] . 
Let us consider the spatially-flat Friedmann-Robertson- Walker (FRW) metric, 

ds 2 = -dt 2 + a{t) 2 (dr 2 + r 2 dn 2 ) (2.2) 



Variation of the action (2.1 ) with respect to the metric leads to the following equations 
(see (201 E3] for details) 



K 2 



H 2 = —(p DE + Pm) (2.3) 
-2H- 3H 2 = k 2 ( PDE + Pm ) (2.4) 
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where 
and 

1 • 



1 rlW 

Pde = -J 2 + V{4>) + Q&FMtf + 24# 3 — ^ (2.5) 

2 dm 



Pde = k 2 



if - V{(j)) - (sH 2 + 2H\ F 1 

- 2H (2FMM + ^) ~ %H 2 ^ - 8H 2 ^ 2 - WHH^i - 16H*f 2 ' 

'(2.6) 

where k 2 = 8ttG. In the present study we assume that the matter sector is modeled by 
an ideal fluid obeying the equation of state (EoS) p m = wp m with constant parameter 
w (mostly non relativistic matter with p m — 0), whose energy density satisfies the 
usual continuity equation p' m + 3Hp m = . The equation of motion for the scalar 
field takes the form 

dV of , „ •■ dF-i -r,\ q 



(j) + m<P + — + 3H 2 [ 2F 1 ((j))(j) + — -<p 2 + im 6 F x {(j))(j)+ 

d<p V d(p J ^ 2 ^ 

YlHHFxUU + 24 (hH 2 + H 4 ) = 

V J dtp 

For generality we keep both Fi(4>) and ^(0) arbitrary, which allows to consider 
different viable cosmological scenarios containing quintessence and phantom phases 
(and early time matter dominated phase, in the presence of matter). As we will 
see bellow, the exponential couplings (that are present in the next to leading order 
corrections to string theory (38], [39], [ID]) appear from reasonable restrictions on the 
relative densities corresponding to kinetic and GB couplings. 



3 The scalar field dominance 



We start studying solutions to Eqs. (2.3 2.6), in the important case when the scalar 



field potential and matter contribution are absent (i.e. V = 0, p m = 0). This leaves 
the two couplings Fi(<p) and ^(0) as the degrees of freedom that may characterize 
the late time cosmological dynamics. Considering the limit of dark energy dominance 
(Pm — >* 0) an d dark energy dominated by the kinetic coupling, we may neglect the 
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free kinetic term and the GB coupling, leading from (2.3) and (2.5) to 



H 2 3F l( f) 2 H 2 



(3.1) 



this limit suggests the following restriction for the kinetic coupling 



F l( j) 2 = k 



(3.2) 



where k is constant. Note that if we divide the Eq. (2.3) by H then with the account 



of (2.5) we may write Eq. (2.3) in terms of the density parameters of DE and matter 
as 

QdE + Q-m — Q<l> + + + = 1 (3-3) 



where 



a 



3H 2 ' 



d F 

tt k = 3F 1 <j) 2 , VL GB = 8H—^<j) = 8H^, Q 



dF 2 p m 
dt ' m ~ 3H 2 



(3.4) 



where we have set k 2 = 1 and p^ = <f) 2 /2 + V, Q k and &gb are the density parameters 
associated with the kinetic and GB couplings respectively. Considering the limit of 



DE dominated by the GB coupling (i.e. — > 0) and , then from (2.3) and (2.5) 
follows 

(3.5) 



which motivates the restriction on the GB coupling 

dFz _g_ 

dt H(t) 



(3.6) 



Then, considering the definitions in (3.4), we may interpret the restrictions on the 



kinetic and GB couplings as the condition that the density parameters fl k and VIqb 
become constants. 

Q k = 3k, Q GB = 8g (3.7) 
where k and g are constants. Having fixed the couplings F\ and F 2 we can solve the 



Eqs. (2.3)-(2.7). 



An important consequence of the couplings of the form (3.2) and (3.6) comes from 
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the generalized Friedmann equation (2.3 2.5). From this equation and making V 
0, follows 

1 



0,Pr. 



H 1 



+ 9kH 2 + 24gH 2 



which imposes the behavior of the kinetic term 



\ 2 H 2 



(3.8) 



(3.9) 



where A 2 = 6 (1 — 3k — 8g), which should be positive, i.e. 1 — 3k — 8g > 0. This also 
indicates that the kinetic term scales in the same way as Pde- Note that if we take 
1 — 3k — 8g = 1, then the important limit <fi 2 = 6H 2 is obtained, which in absence of 
kinetic and GB couplings would correspond to purely kinetic model describing "stiff" 
matter with equation of state w = 1. But the scalar field should satisfy additionally 



the equation of motion (2.7) which depends on the couplings. Multiplying the Eq. 
(2.7) by (f) we can reduce it to a first order equation for the function (p 2 as follows 
(making V = 0) 

1 d(<j) 2 ) 



+ 3H<p 2 + 3H 2 ^(F 1 <p 2 ) + l8H 3 F 1 <p 2 + l2HHF 1 (p 2 + 24:(HH 2 + H 4 

2 at at 



dF 2 
dt 



Replacing F\(j) 2 from (3.2) and dF 2 /dt from (3.6) we find 

1 d{<j) 2 ^ 



2 dt 



+ 3H<j) 2 + 18kH 3 + UkHH + 24gH [H + H 2 ) = 



and taking into account the Eq. (3.9), gives 

(1 - jfe - Ag) H + (3 - 6k - 20g) H 2 



this equation has power-law solution 

// V 



P 



l-k-Ag 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



3 - 6k - 20g 

leading to accelerated expansion provided p > 1. If we take into account the inte- 
gration constant, then this solution could describe phantom behavior with future Big 
Rip singularity. Note that by setting 3 — 6k — 20g = 0, leads to the de Sitter solution 



H = const.. Using (3.13) in Eq. (3.9) we find the scalar field as 

(j) = \p\nt 



(3.14) 
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and the kinetic and GB coupling from (3.2) and (3.6) respectively as 

k 



Ft 



24>/(X P ) 



AV 



A e 2^/(Ap) 

2p 



(3.15) 



where we used (3.14) for t. Note that this exponential behavior is exactly what is 
expected from string theory for the kinetic and GB couplings [38] , [31] (appropriately 



rescaling the scalar field). An important result from Eq. (3.12) is that it can also 
be satisfied when the two constant coefficients are zero. This leads to two linear 
equations for k and g with solution: k = —2,g = 3/4. For this concrete values the 



Eq. (3.11) ((3.12))satisfies automatically for any H(t). This fact can be used to 
reconstruct the model for any given cosmological evolution encoded in H(t). In any 
case if the Hubble parameter differs from power-law, the reconstructed couplings F\ 
and F 2 would be different from simple exponentials. As we have seen the restrictions 



(3.7) are reasonable approximations for late-time cosmology. 



We also can find exact solutions for the equations ((2.3 2.6)) by considering the re- 
strictions 

&k = c k a~ a , Q GB = c g a~ a (3.16) 

where c k and c g are numerical constants. This restriction for the GB coupling has 
been considered in [44J. Introducing the e-folding variable N = In a the equations 



(2.3) and (2.7) take the form 



-H 2 9 + V + 9iJ 4 F^ + 24# 4 ^P 
2 dN 



(3.17) 



and 



1 d 



(H 2 9) + 3H 2 9 
h 12H 



2 dN 

+ 18# 4 (Fi 



dV 

dN 
2 dH 2 dF 2 

~d~N~dN 



J tt2 j 

JFo 



(3.18) 



+ 2AH 4 



dN 







where we have used d/dt = Hd/dN. We have multiplied the Eq. (3.11) by <p and 
represented (d<ft/dN) 2 = 0(N). Using the restrictions (3.16) and additionally consid- 
ering that the scalar field evolves as 9 = const. = X 2 (i.e. = AiV + 0o)> then from 

V = H 2 (3 - A 2 /2 - 3(c k + c g )e- aN ) (3.19) 



((347} follows 



Setting A 2 = 6 to simplify the expression for the potential (i.e. V 

2 ~' !iV ) and replacing in the Eq. (3.18), then the equation for H 2 becomes 

dH 2 



3(c fc + 



c g )H 2 e 



+ 18# z - (c fc 



-c g )e 



-aN 



d\ " v ^ ' 2 

Solving this equation one finds 



dH 2 
~d~N 



+ (2(3 + a)c k + 3(1 + a)c g ) e - aN H 2 = (3.20) 



H* = Ce* N [3c g + 2c k -6e aN ] m 



where C is the integration constant and 

'3c g + 6c fc 



3c„ + 2c k 



+ a 



V2 



24(c g + c fc ) 
(3c g + 2c fc )a 



(3.21) 



(3.22) 



From (3.16) and (3.21) one finds the expressions for the kinetic and GB couplings 



F-i 



18C 



e -( a+Vl )N [ 3Cg + 2ck - 6e aN \ 



-'12 



and 



c g e 



-(a+rn)N 



lC(a + r h )(3c g + 2c k )v2 



2^1 



a + rfi 



6e 



aN 



a 



V2, 



a 3c a + 2cfc 



(3.23) 



(3.24) 



Note that the value rj 2 = leads to the solution known as "stiff" matter H 2 oc e" 



-6N 



Interesting late-time cosmological solutions can be found by setting the power 



772 = 1, or 2. Solving the condition 772 = 1 with respect to a one finds 



H 2 = Ce 



_8(c £ +cj 1 ) 

3c + 2cfc - 6e 3c o+^k 



N 



(3.25) 



Note that all found solutions may be written in terms of the scalar field by replacing 
iV = 0/a/6. In order for this solution to make sense, the integration constant C 
should be negative and 3c g + 2c k < 0. Combining this last inequality with c g + c k < 
(that follows from the positivity of the potential (3.19) for A 2 = 6) leads to the 
restrictions (c g < and c k < —c g ) or (c g > and c k < —3c g /2). Under this 
conditions, the argument of the second exponent in (3.25) becomes negative and the 
sign in the argument of the first exponent depends on the particular values of c g and 



Cfc. Therefore the solution (3.25) is able to explain the late-time cosmic acceleration. 
In terms of the scalar field the potential becomes 

V = l8C(c g + c k )e M - 3C{c g + c k ){3c g + 2c k )e M 



(3.26) 
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where 



= 2 5c g + 2c fc = 2(c g - 2c fc ) 

\/6(3c s + 2c fc ) ' \/6(3c g + 2c fc ) 



Considering 772 = 2, and solving this condition with respect to a gives 



V6cg_ / \/6(e g +e fc ) 
vi*. V I O „ I 0„ „ 3, 1-9,., ^ 



2 



H 2 = Ce ^ a + 2 ^ v \^c g + 2c k - 6e 3c e+ 2c fc v j (3.27) 

where C > and the only restriction on the constants c g and Ck is c g + Ck < 0, to 
keep the positivity of the potential. Therefore the present solution is able to explain 
late-time cosmic acceleration, including super-acceleration (w e // < — 1) if one of the 



arguments in the exponentials in (3.27) is positive. The scalar potential takes the 
form 

V = Voe^ (3c g + 2c k - 6e 72 *) 2 (3.28) 
where V = -3C(c g + c k ), 7l = 72 = -^g^. 

Stability of solutions 



To board the issue of stability of the solutions corresponding to the restrictions (3.7) 



and (3.16) we will introduce the autonomous system for the Eqs. (2.3)-(2.7) in absence 
of matter. We used this approach because of the nature of the restrictions we have 
considered, that have direct connection with the dynamical variables defined for this 
model. The stability analysis based on perturbations, for inflationary solutions in 
more general second-order scalar tensor theory, has been performed in [33]. Let's 
consider the following dynamical variables 

K,(J) 

2 • at 2 ti 
3 = 8kH ^ e = lP 
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The terms with derivatives in (2.3)-(2.7) transform under the change of variable as 



V2H 2 {xe + x') 



K 



dV H 2 (2ye + y') 



d 2 F 2 



\/2kx 
xg' — g{2xe + x') 
~d~4? ~ lQH 2 x 3 
dFi xk' — K>k(xe + x') 



(3.30) 



3V2H 



2 X A 



where "prime" denotes the derivative w.r.t the e-folding variable N. In the subsequent 
analysis we will set k 2 = 1. In terms of the variables (3.29) and using (3.30), the Eqs. 



(2.3)- (2. 7) can be transformed into the following first-order autonomous system 



x 2 + y + 3k + 3g -3 = = 
2xx' + 2(3 + e)x 2 + y' + 2ey + k 1 + 2(3 + 2e)k + 3(1 + e)g = 



(3.31; 



2e + 3 + x 2 - y 



-k' 



; (3 + 2e)k- g' - (2 + e)g = 



We should consider that the solution (3.13 )-(3.15 ) was obtained without potential (i.e. 



V = in (3.29) and therefore there is not variable y associated with the potential. 



Additionally we have to take into account the constraints on the dynamical variables 



k and g (3.7), that lead to the solution (3. 13)- (3. 15). All this amounts to complement 



the dynamical system (3.31) with the following equations 



y = 0, g' = 0, k' = 



(3.32) 



And therefore the system reduces to the following one- dimensional dynamical system 
for x 

3g + 3k + x 2 -3 = (3.33) 
3g(e + 1) + 2e{2k + x 2 ) + Qk + 2x(3x + x') = (3.34) 

where 



3 (2g + k - x 2 - 3) 
3g + 2k-6 



(3.35) 
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Note that the constants k and g are not arbitrary, as they should satisfy the restriction 



(3.33). From (3.34) and (3.35) follows 



x 



3(g-2x 2 )(3(g + k-l)+x 2 
2x{3g + 2k - 6) 



(3.36) 



The critical points of (3.36) are: x\^ = ±>Jg/2 and x^^ = ±a/3(1 — g — k). It is 
easy to check that the solution (3.13), (3.14) is a critical point provided A = \/2xi^ 



or A = -\/2^3,4- Taking the derivative of (3.36) w.r.t. x and evaluating at the critical 
points, we find the eigenvalue equations 



1. Ai 



dx I x — x l,2 



- 3~^j^Z6 • The stability of the fixed point x\p demands that 
Ai < 0, which leads to one of the following conditions: (g > 6, k > |(6 — 7g)), 
(g < 6, k > |(6 - 3(7)), (g < 6, fc < |(6 - 7#)) or (<? > 6, k < |(6 - 3^)). If one 
of this conditions is satisfied, then the fixed point x\^ is an attractor. 



2. A, 



4f I 

dx Iz=:r3,3 



3g+^^6 • ^ or -^2 < this fixed point is an attractor, which 



leads to the conditions: (g < 6, |(6 — 7g) < k < |(6 — 3g)) or (g > 6, 
|(6 -3p) <fc< |(6-7^)). 

Taking into account that the solution is a fixed point x = xq = A/v^2, then A takes 
the following values: A = ±y^ for the fixed points x\^ and A = ±a/6(1 — k — g) for 



X34. Using this results in the Eq. (3.33), then the above inequalities translate into 



the following conditions on g to achieve the stability of the solution (3.13 )-( 3.15 ): 



Ai < : < g < 12 



A 2 < : g < or g > 12 



(3.37) 
(3.38) 



For the second solution we have used the restrictions (3.16) on the density parameters 



for the kinetic and Gauss-Bonnet couplings as defined in (3.4). Note that these density 



parameters coincide with the ones defined as the corresponding dynamical variables 



in (3.29) (i.e. k = and g = Q g , setting k 2 = 1). In addition, to obtain the solution 



(3.21) we have introduced the ansatz: 



s'2 



(3.39) 
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where we have limited to the case A 2 = 6, which gives the potential of the form 



V = V H 2 e 



-aN 



V = -3(c fc + Cg 



(3.40) 



Form Eq. (3.39) follows that 



<P = 0o + VQN (3.41) 

And therefore the variable x becomes constant as follows from the definition x = 
-J^ = 4| = \/3 = Xq. This leads to x' = 0. Setting O — 0, we can re express 
the variable N in terms of <f), and the potential takes the form: V = VoH 2 e~ a< ^^, 



which according to (3.29) gives y = V^e "^/^ and from (3.16): k = c^e "^/v^ 



g = c g e ^l^. Using this results, the dynamical system reduces to the equations 



a 



a 



y 



k' 



a 



xk 



(3.42) 



From these equations follows that the only critical point is (0, 0, 0) and the eigenvalue 
matrix becomes diagonal with all the eigenvalues equal to 

a 



Ai — As 



A, 



-a 



(3.43) 



from which follows that the point (0, 0, 0) is a stable critical point (attractor) if a > 0. 
This makes sense if we see that at t — > oo (or a — > oo) (as follows from the solution 



(3.21) for the discussed restrictions on Cfc and c g ), then g — >■ 0, k — > and y — > for 



positive a, as follows from (3.16) and (3.40). 



Stability under tensor and scalar perturbations 

The present model is a particular case of the generalized Galileon theory 
|47j . [48J. The conditions to avoid ghost and gradient instabilities in the generalized 
Galileon theory have been presented in [45]. Here we use the stability analysis pre- 



sented in (45] applied to the model (2.1) as special case in absence of matter. The 
generalized Galileon model may be written as [45J (setting k 2 = 1) 



S 



d x\J—g 



K{<j>, X) - G 3 (0, X)D<f> + G 4 (0, X)R + ((D0) 2 - (V M V,0) 2 ) 



+ G 5 (<f>,X)G flv V}j,V u <f> 



l dG 5 
6 dX 



[(□0) 3 - 3(D0)(V At V,0) 2 + 2(V / ,V,0) 3 ] 



(3.44) 
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where X = -V M 0V>/2 and (V^V, 



(V M V^0)(V^V CT 0)(V (T V^0) . To obtain 



the model (2.1) (with S m = 0) we use the following correspondence: 



K(6,X) = -V(<f>)+X- 



d A F 



2 V 2 



d 3 F 2 



X 2 (3-lnX), G 3 = -A— fX(7-31nX) 



Ga 



1 .d 2 F 2 



-4— ^X(2-lnX), G 5 



+ 4^1nX 



(3.45) 



1 + 



where the first term in the expression for G 5 gives (up to total derivative) the non- 



minimal kinetic coupling as appears in the third term in (2.1) (note that Fx — 1 leads 
to G5 = —(f)), and all the terms depending on F 2 in (3.45) reproduce the GB coupling 



in (2.1). 



Tensor perturbations 

The quadratic action for tensor perturbations hij is given by (see [43] ) 

1 



5 2 Sn 



dtd 3 xa 3 ( Q T h\ - ^f(V^) 2 
a 



(3.46) 



where 



r 2x dGi x(h^ g * dG " 



Ga — X 



■dG 5 dG 5 
' dX 86 



(3.47) 



from (3.46) follows that the conditions to avoid ghost and gradient instabilities under 



tensor perturbations reduce to 



Gt>0, J r T >0 



(3.48) 



Replacing the solutions (3.13 )-(3.15 ) in (3.47), taking into account the definitions 



(|3.45|), we can analyze the conditions (|3.48|). The condition Q T > takes the from 

(3.49) 



fc(AV + 2A P + 40 2 ) ^ 

SJt — 1 — °9 ^r~r^ > 



(Xp + 20) 2 

Then the tensor ghost is absent if the above inequality is satisfied, or if g < and 
k < (provided 6 = Aplnt > 0). Note that at large times Qt behaves as 



Sg-k 



(3.50) 
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which is also the de Sitter limit that takes place at p — >■ oo (in this limit the EoS 
parameter w = — 1 + 2/(3p) takes the value w — — 1). So the condition k < 1 — 8g 
guarantees the absence of ghosts instabilities at large times and in the de Sitter limit. 
The condition for Tt is given by 



1 8o fc(AV + 2A P + 40 2 ) 
:Ft - 1 ~J + (Xp + 2^ >0 

at large times and in the de Sitter limit Tt behaves respectively as 

1 % 



T 



P 



1 + jfc 



(3.51) 



(3.52) 



Combining with the previous conditions, the avoidance of ghost and gradient insta- 
bilities under tensor perturbations impose the following restrictions at large times: 



P > 0, 9 < 



p 



ig-p 



< k < 1 — 8g 



4 + 4p p 

and at the de Sitter limit we find the following restrictions 

1 



9 < 



-1 < k < 1 - 8g 



(3.53) 



(3.54) 



Anther interesting result is obtained in the case when we neglect the kinetic coupling. 



By setting k = in (3.49) and (3.51) it follows that Qt and Tt become constants 
equal to 

8q 

Q T = 1 - 8g, T T = 1-— (3.55) 



Thus in the particular case without kinetic coupling, the model is free of ghost and gra- 
dient instabilities (independently of time) under tensor perturbations, for the power- 
law evolution, provided g satisfies the restrictions 



0<p<l, g<p/8, or, p > 1, g<l/8 



(3.56) 



Let's consider the second solution given by Eqs. (3.19) (with A = 6), (3.21 )-(3.24) 



in the N variable. To evaluate the conditions (3.48) for this solution, we use d/dt 



Hd/dN in Eqs. (3.45) and (3.47). In this case the expressions for Qt and Tt are 
more involved, but assuming a > we have found that at large times (future infinity) 
Qt — ^ 1 an d Tt — > 1 (independently of c g and c^), which indicates that we can 
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expect stability under tensor perturbations at future times. On the other hand, if we 
consider a < then the stability conditions take place at early times. Nevertheless, 
these appreciations correspond to asymptotic behavior. To give a numerical example 



we analyze the specific viable model of DE (3.27), (3.28) corresponding to r/2 = 2. To 



simplify the numerical analysis we use the units k 2 = 
we will assume that the current value of the EoS is w 



1,Hq = 1. As initial condition 



-1 and the constants c g , Ck 
and C are subject to the flatness condition (here the constant C is given in units of 
Hq). Additionally we can restrict the constants c g and Ck in such way that the "stiff" 
matter term (cx a -6 ) in the expression for H disappears. This gives the relation 



3(g 2 -8g) 
2(9-6) • 



The numerical analysis for the specific values of c q = 10 and 



Ck = —2 x 10 12 is resumed in Fig.l. For these values the DE model describes with 
high accuracy a de Sitter universe, with EoS practically constant, w = —1. 




Ln(l+z) 



Fig. 1 Qt and Tt versus ln(l + z), for c g 



10 



-12 



and Ck = — 2 x 10 12 . Note that 



Qt and Tt are practically constants equal to 1, during all cosmological epochs. The 

internal graphic shows the future behavior of Qt and Tt for the redshift z in the 
region (—0.99, —0.94). Note that starting form z ~ —0.97 the values of Qt, Tt begin 
to differ from 1. This solution is stable under tensor perturbations for a wide range 

of redshifts 

So we see that the solution is stable under tensor perturbations during different cosmo- 
logical epochs, including late times and future evolution, and starting at z ~ —0.98, 
Tt tends to negative values, signaling instabilities at far future when the Eos tends 
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to w\ z ^.-i = —3. 
Scalar perturbations 

A similar expression for the quadratic action for scalar perturbations ( is given by 



S 2 S S = J dttfxa 3 (W 2 - J(VC) 2 ) (3.57) 



where 



(3.58) 



with E and 6 given by 



dX ' dX 2 ' V" 5X ' " 5X2 7 V d <t> ' "<90<9X") 



(3.59) 



(3.60) 

The ghost and gradient instabilities are absent provided 

Qs > 0, > (3.61) 
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In the case of the solution (3.13 )-( 3. 15 ), and using (3.45) and (3.58)-(3.60), we find 
g s = 3 (A(8g + k - 1) ln 2 (t) + 2(16g + k - 2) ln(t) 

+ 8g + k - l) (l6(64# 2 + 8g(4k - 1) - Qk 2 + 3Jfe - 1) ln 4 (t) + 8(256# 2 + 8#(9£; - 4) 
- 15A; 2 + 12k - 4) ln 3 (t) + 4(384/ + 80gk - 48g - 21k 2 + 18k - 6) ln 2 (t) 
+ 2(256# 2 + 8g(9k - 4) - 15A; 2 + 12k - 4) ln(t) + 64g 2 + 8g(4k - 1) - 6k 2 + 3Jfe - l) 
/ [(2 ln(t) + !) 2 {4{3k - 1) ln 2 (t) + (6k - 4) ln(t) + 3k - l) 2 ] 



and 



3{-Ug 2 {p + 1) + 8g{-2kp + k + 2p + 1) - k(4kp + k - 2) - 1) 



k{4p + 3) 



3p(3k - 1) 

8(1 - 12g) 2 k((p + l)ln(t) 



(3.62) 

4kp + k 
3p(21n(t) + l) 



+ 



3p(2 ln(t) + l) 2 3p(3A: - 1) [2 ln(t)((6fc - 2) ln(t) + 3k - 2) + 3k - 1} 
4k (16(1 - 12g) 2 k{(3k - 2) ln(t) + 3k - 1) 



3p(21n(t) + l) 3 3p{3k - l)[21n(t)((6fc - 2) ln(t) + 3A; - 2) + 3A; - l] 5 



At large times, taking the limit t — > oo in (3.62) and (3.63) we find 

3(8g + k — !){<o4g 2 + 8g(4k - 1) - 6k 2 + 3k - 1) 



IS t-KX3 



and 



(1 -3fc) 2 

1 + 2Jfe - 64# 2 (1 + p) - A; 2 (l + 4p) + 8g{l + k + 2p - 2kp) 



(3.63) 
(3.64) 

(3.65) 



(-1 + 3% 

So the stability under scalar perturbations at large times is possible, provided Gs\t-+oo > 



and J-s\t-+oo > 0. Note that the same limits as described by Eqs. (3.64) and (3.65) 
take place at early times t — y 0, and therefore the same conditions guarantee the 
stability at early and late times. 



Taking the limit p — > oo in (3.65), one finds the de Sitter limit 

4g + 4gk + k 2 ) 



S\dS 



4(lQg 2 



(3.66) 



1 - 3k 

Note that Gs\t^oo does not depend on p. Therefore the solution is stable under scalar 
perturbations in the de Sitter limit if g and k satisfy the restrictions 

(8g + k-l)(<64g 2 + 8g(4k-l)-Qk 2 + 3k-l) > 0, 16# 2 -4c/+4c/A; + A; 2 > 0, l-3k > 

(3.67) 
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If we neglect the kinetic coupling (k = 0), then from (3.62) and (3.63) one finds that 
Qs and Ts are constants given by 

1 + 64# 2 (1 + p) - 8#(1 + 2p) 



Qs = 3 - 38V + 1536<? d , F s 



V 



(3.68) 



Thus, in this case the model is free of ghost and gradient instabilities (independently 
of time) under scalar perturbations, for the power-law evolution, provided g satisfies 
the restrictions for p > 



3 - 384g 2 + 1536# 3 > 0, 1 + 64# 2 (1 + p) - 8g(l + 2p) > 



(3.69) 



The accelerated expansion takes place for p > 1, which together with (3.69) leads to 
the restrictions: 
For 1 < p < 1 + y/5 



1 , t . 1 + 2p 1 / 4» 2 - 3 1 . t . , 

i6 (1 -^» <s< i6(T^-T6V(TT^' or ' »>^ 1 + ^ (X70) 



For p > 1 + y/h 
(l-VE)<g< 



l + 2p 1 Up 2 -3 



16(1+ p) 16\/(l+p) 2 ' 



or, # > 



1 + 2p 1 Up 2 -3 



16(1+ p) 16\/(l+p) 2 



(3.71) 



Turning to the second solution (3.21), we considered the case rj2 = 2, but the ex- 
pressions for Qs and J-g are long enough to be displayed here. We performed the 
numerical analysis for the same particular case considered for the tensor perturba- 
tions. In Fig. 2 we resume the numerical results obtained for the study of stability 
under scalar perturbations. 
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Fig. 2 Qs and JFs versus ln(l + z), for c g = 10~ 12 and Ck = —2 x 1CT 12 . It was 
found that Q s = 3 and J^s = with high accuracy during a wide range of redshifts 
which covers the epochs relevant for dark energy. The internal graphic shows the 
future behavior of Qt and Tt for the redshift z in the region (—0.99, —0.97). Note 
that starting form z ~ —0.99 the values of Qs, Ts begin to change. 

Numerical results show that after z ~ —0.998 the quantity Qs becomes negative and 
Ts starts increasing to large positive values. In summary, we have found that the 
solution is free of ghost and gradient instabilities under scalar perturbations over a 
wide range of redshifts, covering different cosmological epochs, including those that 
are relevant for DE. 



4 Scalar field and matter 



Let's consider the complete model with the contribution of matter sector. Combining 
the equation of state p m = wp m with the continuity equation p m + 3H(p m +p m ) = 0, 
leads to the energy density of matter p m = p m o e_3( ' 1+w ^- It is usually assumed that 
the matter sector is composed mostly of non relativistic matter with w = 0. In the 



e-folding variable N, the equation (2.3) takes the form 

,dF 2 



H 1 



-H 2 9 + V + 9H 4 Fx6 + 2AH 4 — 
2 dN 



+ Pmoe 



-3(l+w)N 



(4.1) 



and the equation of motion (2.7) remains as Eq. (3.18). Defining the same limits as 



before to determine the kinetic and GB couplings, in the e-folding variable take the 
form 

k dF 2 g 



(4.2) 



dN H 2 

Using the definition of the density parameters, this corresponds to cosmologically 
viable limits for the density parameters = const. = 3k and VLqb — const. = 8g. 



It remains to determine the potential in order to integrate the equations (4.1) and 
(3.18). If we demand that Qfp = const., then Qde = const, and independently of 



the potential the solutions correspond to scaling regime where the energy density of 
the DE mimics the matter density p m . In fact the condition fL = const, can be 
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achieved by making <fi 2 oc H 2 and V oc H 2 . But to have a transition to accelerated 
expansion, we may opt for restrict the potential in such a way that the kinetic term 
in the Friedmann equation (4.1) behaves like (f) 2 oc H 2 . The restriction 

V + p m0 e- 3(1+ - )Ar 



3aH 2 



(4.3) 



reduces the Eq. (4.1) to the form 

9 = 6 (1 - a - 3k - 8g) = 9 (4.4) 
where 1 — a — 3k — 8g > 0. Integrating this equation gives the scalar field as <fi = 



\/9qN + 0o- The presence of matter term in (4.3) marks the difference with scaling 
behavior. The limit V — > gives rise to scaling solutions, and we should expect 
that at early times the scalar potential is negligible or is proportional to p m . If we 
introduce the relative potential density Qy — V/(3H 2 ), then the restriction (4.3) is 
equivalent to Qy + Q m = a. Conversely, we may consider (4.4) as an ansatz for 
the free kinetic term (i.e <ft 2 /(6H 2 ) = const.) (see [44J) and arrive at the form of the 



potential (4.3) from the Friedmann equation (4.1). As will be seen the restrictions 



(4.2) and (4.3) (or (4.2) and (4.4)) are good approximations for late-time cosmology. 



Replacing (|4.2|)-(|4.4|) in the equation of motion (3.18) one obtains 

3a -Qk- 20g) H 2 + (1 + w)p m0 e^ 1+w ^ N -- 
k — 4g = or 3 — 3a — 6k 



jtt2 

'l-k-Ag) — + 2(3 



(4.5) 
- 20g = 0, 



Satisfaction of one of the two conditions 1 
leads to the solution H 2 oc e~ 3 ( 1+w ^ N which according to (3.3) imply scaling behavior 



of the DE. Solving (4.5) gives 

H 2 = Hi { V n m0 e-^ N + n oe -r") 



(4.6) 



n 



7 



3a + 6fc + 20#-3 



6a + 9£; + 2% + 3u;(l-£;-4#)-3' ' k + Ag - 1 

where Qo is the integration constant and Q m o = p m o/(3-ffo) (-^o is the current value 
of the Hubble parameter). Note that the term nQ m Q can be interpreted as an effective 



matter parameter. The restriction (3.3) imply that rjVt m Q + fi 
scalar potential takes the form 



1. From (4.3) the 



V = (ar, - l) Pm0 e- 3(1+w)7V + 3aH 2 Q e^ N 



(4.7) 
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Replacing the solution (4.6) into Eqs. (4.2), one finds the expressions for the kinetic 
and GB couplings as 

k 



(4.t 



where 9 is given by (4.4). 



7 



3w - 2 7 + 3 rjVt m0 e-^ w ~^ N 



3w — 7 + 3 3w — 7 + 3 



fin 



(4.9) 



An interesting case takes place in Eq. (4.7) when ar) = 1. In this case the matter 



term in the potential disappears, giving place to the simple exponential behavior 

V = 3aHlVL Q e~ lN = 3a^f2 e-^ /V5 ° (4.10) 

where —3a = 9k + 28g + 3u>(l — k — 4g) — 3, and we used the solution for the scalar 
field with 0o = 0. According to the solution (4.6) and Friedmann equation (2.3) the 
DE density is given by 

Pde = (v~ l)p m0 e- 3{1+w)N + 3n H 2 e-' N = { V - l)p m0 (l + z) 3 ^ + 3fi # 2 (l + z)f 

(4.11) 

where we used the redshift relation: N = — ln(l + z). It is clear that the scalar field 
keeps track of the DM at early times if 7 < 3(1 + w) (or more sharply if 7 < 0, 
since in this case the second term becomes negligible for z » 1). This solution 
agree with the fact that during the matter dominated epoch wde ~ w. The first 



term in (4.6) is an adjustment of the initial parameter appearing in the Lagrangian 



density of the matter due to the interacting terms, which we interpret as an effective 
matter parameter. Therefore, is the quantity r]Q m0 which has to be adjusted with the 
observational data and should be about 0.27. If rj = 1, the DE presents power-law 
behavior with constant EoS wde = — 1 + 7/3. The condition 7 < in (4.11 ) leads to 
quintom behavior, since it allows the transition to the phantom phase. 
In the important case of non relativistic dark matter component with w — 0, the DE 



equation of state that follows from (4.11) is 



w DE 



[1 - r/^ m0 )(3 - 7)e' 



-7JV 



3(1 - rAno)e-^ + 3(77 - l)Q m0 e- 



3N 



(4.12) 
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the current value wde(0) = wq may be written as 

1 - rjVt m0 



w 



1-Q 



771.0 



1 + 



7 



(4.13) 



The first factor in this equation qualifies the difference with respect to constant EoS, 
which is the case for 77 = 1. According to current observations wde could be very 



close to — 1; so setting w = — 1 in (4.13) we obtain the relation 



7 



3(1 - r])fl mQ 



(4.14) 



1 - i]Q mQ 

On the other hand, the deceleration-acceleration transition occurs at redshift zt when 
w e ff = —1/3. Evaluating w e ff = — 1 — |^ from (4.6), and setting w e ff(zr) = —1/3 
we find another relationship 

3(2-7) 



V 



3(2- 7 )+ 7 (l + z T ; 



3-7 



(4.15) 



where we used (4.14). Thus for example, taking = 0.7 and 7 = —0.1, gives r\ w 1.09 



and from (4.14), f2 m0 f» 0.26. So, the restrictions (4.2) and (4.3) allow to describe a 



viable cosmological scenario. 

We can make another interesting choice for the kinetic and GB couplings, namely we 
consider the ansatzes for the density parameters 



c k e 



-XN 



n 



GB = c„e 



-XN 



(4.16) 



this means that f^, Qgb scale as a~ x , where A should take such value that guarantees 
the early time dominance of the matter sector and lower rate of decaying of dark 
energy during expansion at late times. With this choice the Friedmann equation 



(4.1 ) becomes 



3H 2 



-H 2 e + V + 3(c fc + c g )H 2 e- XN + p mQ e 



(4.17) 



and restricting the potential in the form 



V + 3(c k + c g )H 2 e~ XN + Pm0 e- 3 ^ N = 3aH 2 , 



(4.18) 
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reduces the Eq. (4.17) to: 9 = 9 = 6(1 — a), where < a < 1. Replacing (4.16), 



(4.18) and 6 in the equation of motion (3.18) one obtains 

jtt2 

(3 - (c k + 3c g /2)e~ XN ) — + (18(1 - a) + [2(3 + X)c k + 3(1 + X)c g ] e~ XN ) H 2 



+ 3(1 + w)p m0 e 



-3{l+w)N 







(4.19) 



In order to simplify the integration of this equation we take a 
and considering the important case of non relativistic presureless matter with w 
we find 



24(c g +c fc )+3c g A+2e fc A 
6(3c s +2c fe ) ' 

= 0, 



H 2 



6a lPm0 (l + z) 3 , C(l + z)- 2 ^+ x 



+ 



(6 - ax{\ + z) x )a 2 6 - a x {l + z 
where we used the relation: N = — ln(l + z), and a± = 3c g + 



(4.20) 



2c fc , a 2 = 3c 9 (A + 5) + 



2cfc(A + 9) and a 3 = 3c 9 (A + 1) + 2c^(A + 3). Using this result and Eqs. (4.2) we find 
the following expressions for the kinetic and GB couplings 



where #0 = 6(1 — a 
2 (3-A)a 2 



F x 
), and 
2F1 



6 — (Tie 



— XN 



601 

C"2 



p m0 e- 3iV + Ce( 2CT 3M-A)JV : 



(4.21) 



A-3 3 Qe XN 
A A o"i 



2^ 



'^1 



1 1 4- 2 °" 3 66 



Ao"i 



(4.22) 



Note that the solution (4.20) has sense only if A < 0. In this case at early times (z — > 
00), the Hubble parameter behaves as H 2 ~ (a 1 /a 2 )pm (l + z) 3 + (C/6)(l + z) x - 2 ' 7 ^' 72 . 
The first term is the usual matter term and the second term behaves as (l + z) x ~ 2aa / a2 . 
So to maintain matter dominance at early times is sufficient the condition \—2a^ja2 < 
2, which satisfies the nucleosintesis constraint (with appropriate constant C) on dark 



energy in early cosmology. An interesting fact of the solution (4.20) is the existence 
of future singularity at the redshift z s = (6/o"i) 1 / A — 1, where o\ should satisfy the 
inequality < (6/cri) 1 / A < 1 in order to have the singularity in the future. Is easy to 
see from (4.20) that at this point H — > 00, which together with the singularity in the 
Hubble parameter signify that the energy density and pressure become infinity at z s . 
By definition, at this point a — l/(z s + 1) is finite, so we have a type III singularity 
J). Note also that from Eqs. ( |4.21| ) and ( |4.22[ ), it follows that the kinetic 
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coupling disappears (Fx —> 0) and the GB coupling diverges at the singular point z s . 



From Eq. (4.20), the DE equation of state is given by (assuming presureless matter) 

WDE = Jl Ljk 42 3 

3H 2 -3Q m0 (l + z) 3 

where H = H/Hq and Hq is the current value of the Hubble parameter. As an 
example, a set of parameters satisfying all discussed above restrictions is: U\ = 4, cr 2 = 
18, cr 3 = 1, A = —1/4, Q m0 = 0.27, C = 0.92, and the constant C is related to the other 
parameters through the flatness condition: C = C/Hq = 6 — <j\ — 18aiQ m o/a2- For 
this set of parameters, at early times wde\z^oo = 0, which means that at early times 
the DE mimics the usual matter behavior (or early time matter dominance as follows 



from (4.20). The current value of DE EoS is u>£>-e(0) ~ —1.02, and the asymptotic 
value at far future is wde\z->oo = —1 — 1/6. The type III future singularity takes place 
at z s ~ —0.8. Note that although the type III singularities are weaker than big rip 
singularities, and in some cases one can extend the evolution behind this singularity 
[5Tj . [52"] . [53], in the concrete case considered here, although the density becomes 
finite after the singularity, it also becomes negative (i.e H 2 < for t > t s ), which 
makes it impossible for the universe to evolve beyond this singularity. Therefore, one 
can distinguish an scenario of the singularity of type III, which is stronger in the sense 
that the evolution beyond it, is not physically allowable. 



5 Discussion 

We studied late time cosmological solutions based on string spired scalar-tensor model 
including a coupling to the Gauss-Bonnet invariant and kinetic couplings to curva- 
ture. We presented various exact cosmological solutions that exhibit the deceleration- 
acceleration phase and even the very probably phantom phase. These phases are 
necessary for the successful explanation of the early inflation as well as the currently 
observed phase of acceleration (super-acceleration) in late universe. 
As a criteria to fix the couplings we have based on asymptotic limits of the model. 
Thus in absence of potential, the limit of vanishing time derivative of the scalar field 
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(0 ~ 0) suggests the constraint (3.6) imposed on the GB coupling. Neglecting the 
kinetic term, the limit (F 2 — >• 0, V — > 0) suggests a cosmologically viable restriction 



on the kinetic coupling (3.2). The limit (Fi — > 0, F 2 — > 0, V — > 0) suggests a behavior 



of the kinetic term (3.9) which is sufficiently general to obtain interesting cosmolog- 
ical scenarios. The restrictions on the interaction terms are also equivalent to the 
constancy of the respective density parameters, i.e. = const, and Vt GB = const. 



(see (3.7)), that are reasonable limits for a late time universe. These restriction in the 
scalar field dominance case, lead to the exponential couplings characteristic of the low 
energy a' expansion in string theory [38], [39], [3D]. In the more realistic case with 



the presence of matter, the restrictions (3.7) (or (4.2)) lead to solution (4.6) with 
dark energy that keeps track of the earlier matter dominance, and containing late 
time quintessence and phantom phases. In this case, to maintain the behavior of the 
kinetic term as oc H 2 , we restrict the scalar potential according to the expression 



(4.3), which leads to the exponential behavior (4.7, 4.10). 

Another cosmologically viable expressions for the kinetic and GB couplings are ob- 
tained by imposing scaling behavior of the type a~ x to the corresponding density 
parameters, i.e. Qk oc a ~ X an d Qde oc a~ x , where we considered the same power in 
order find exact analytical solutions. In the case of scalar field dominance, we found 



a cosmo logical solution (3.21) able to explain late time accelerated expansion and 
even a transition to phantom phase. We considered the autonomous system for the 
model in absence of matter and have found the conditions for stability of the solutions 



(3.13) and (3.21). We have also studied the stability properties of the solutions (3.13) 



and (3.21) under tensor and scalar perturbations of the metric. For the power-law 
solution we have found the restrictions on the parameters in order to avoid ghost 
and gradient instabilities in the following cases: at large times, in the de Sitter limit, 



and in absence of the kinetic coupling. For the solution (3.21) we considered a DE 
scenario very close (with high accuracy) to the de Sitter solution, and it was proved 
that the solution is free of ghost and gradient instabilities under tensor and scalar 
perturbations, from early times to the present and into the future up to z ~ —0.99. 
In the model with matter content, a variety of exact solutions can be found depending 



on the correlation between the parameters of the equation (4.19). In fact, the general 



26 



solution to the cosmological equation (4.19) is given through hyper-geometric func- 



tion, but using the limits when this function simplifies to elementary exponentials we 
found relations between the parameters. For one particular choice of the parameters, 



we found the solution (4.20), which unifies early time matter dominance with late 



time accelerated expansion including transition to phantom phase. This solution also 



presents future type III singularity provided <T\ in (4.20) satisfies < (Q/c>i) l l x < 1. 
An specific characteristic of the singularity that appears in the solution (4.20) is that 
the density becomes negative for times t > t s making it impossible for the universe 
to evolve beyond this singularity. We cited one example, but there are another re- 



lations between the parameters of Eq. (4.19) that could realize phenomenologically 



acceptable cosmological scenarios. In all considered solutions the scaling behavior of 
the kinetic energy is the same as of the total energy density, i.e. 2 oc H 2 , which 
in absence of interactions corresponds to a conformal limit of the model (a purely 
free kinetic scalar model is conformal invariant). We have studied the dark energy 
problem in the context of the scalar-tensor model, with kinetic and GB couplings con- 
tributing to the total energy density. By imposing physically motivated restrictions 
on density parameters involving the interaction terms, we significantly simplified the 
cosmological equations and have obtained exact viable solutions, compatible with the 
current phenomenology of the dark energy. 

References 

[1] A.G. Riess, et al., Astron. J. 116, 1009 (1998); astron. J. 117, 707 (1999). 
[2] S.Perlmutter et al, Nature 391, 51 (1998) 

[3] M. Kowalski, et. al, Astrophys. Journal, 686, p.749 (2008), arXiv:0804.4142 

[4] M. Hicken et al, Astrophys. J. 700, 1097 (2009) [arXiv:0901.4804 [astro-ph.CO]]. 

[5] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 
arXiv:0803.0547 [astro-ph]. 



27 



[6] W. J. Percival et al., Mon. Not. Roy. Astron. Soc. 401, 2148 (2010); 
arXiv:0907.1660 [astro-ph.CO] 

[7] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003) [arXiv:astro- 
ph/0207347] 

[8] T. Padmanabhan, Phys. Rept. 380, 235 (2003) [arXiv:hep-th/0212290] 

[9] V. Sahni, A. Starobinsky, Int. J. Mod. Phys. D9, 373 (2000); arXiv:astro- 
ph/9904398. 

[10] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D15 1753 (2006), 
arXiv:hep-th/0603057 

[11] V. Sahni, Lect. Notes Phys. 653, 141-180 (2004), arXiv:astro-ph/0403324v3 

[12] T. Padmanabhan, Phys. Rept 380, 235 (2003), [hep-th/02 12290]. 

[13] B. Boisseau, G. Esposito-Farese, D. Polarski, A. A. Starobinsky, Phys. Rev. Lett. 
85 (2000) 2236; arXiv:gr-qc/0001066 

[14] L. Perivolaropoulos, JCAP 0510, 001 (2005); arXiv:astro-ph/0504582 

[15] Fujii Y and Maeda K I 2007 The Scalar- Tensor Theory of Gravitation (Cam- 
bridge Monographs on Mathematical Physics) (Cambridge University Press) 
ISBN 0521037522 

[16] S. Nojiri, S. D. Odintsov, Phys. Rept. 505, 59 (2011); arXiv:1011.0544[gr-qc] 

[17] L. N. Granda, JCAP 07, 006 (2010); arXiv:0911.3702 [hep-th] 

[18] L. N. Granda and W. Cardona, JCAP 07, 021 (2010); arXiv: 1005.2716 [hep-th] 

[19] S.V. Sushkov, Phys. Rev. D80, 103505 (2009); arXiv:0910.0980 

[20] L. N. Granda, Class. Quant. Grav. 28, 025006 (2011); arXiv: 1009.3964 [hep-th] 

[21] E. N. Saridakis, S. V. Sushkov; Phys. Rev. D81, 083510 (2010); 
arXiv:1002.3478v2 [gr-qc]. 

28 



[22] C. Gao, JCAP 06, 023 (2010); arXiv:1002.4035[gr-qc]. 

[23] L. N. Granda, E. Torrente-Lujan, J.J. Fernandez-Melgarejo, Eur. Phys. J. C 71, 
1704 (2011); arXiv: 1106.5482 [hep-th] 

[24] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D71, 123509 (2005); hep- 
th/0504052. 

[25] T. Koivisto and D. F. Mota, Phys. Rev. D75, 023518 (2007); hep-th/0609155 

[26] T. Koivisto and D. F. Mota, Phys. Lett. B644, 104 (2007); astro-ph/0606078 

[27] S. Tsujikawa and M. Sami, JCAP 0701, 006 (2007) [arXiv:hep-th/0608178] 

[28] B. M. Leith and I. P. Neupane, JCAP 0705, 019 (2007) [arXiv:hep-th/0702002] 

[29] R. A. Brown, R. Maartens, E. Papantonopoulos and V. Zamarias, JCAP 0511, 
008 (2005); [arXiv:gr-qc/0508116] 

[30] S. Nojiri, S. D. Odintsov, Phys. Lett. B631, 1 (2005); hep-th/0508049 

[31] S. Nojiri, S. D. Odintsov, M. Sami, Phys. Rev. D74, 046004 (2006); hep- 
th/0605039 

[32] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, S. Zerbini, Phys. Rev. D73, 
084007 (2006); hep-th/0601008 

[33] S. Nojiri, S. D. Odintsov, and S. Ogushi, Int. J. Mod. Phy. A 17, 4809 (2002). 

[34] I. P. Neupane, B. M. N. Carter, JCAP 0606, 004 (2006); hep-th/0512262 

[35] S. Nojiri, S. D. Odintsov and P. V. Tretyakov, Phys. Lett. B651, 224 (2007); 
arXiv: 0704. 2520 [hep-th] 

[36] S. Nojiri, S.D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007) 

[37] L. Amendola, C. Charmousis and S. C. Davis, JCAP 0612, 020 (2006); hep- 
th/0506137 



29 



[38] R. Metsaev and A. Tseytlin, Nucl. Phys. B 293, 385 (1987). 

[39] C. Cartier, E. J. Copeland, R. Madden, JHEP 0001, (2000) 035; hep-th/9910169 

[40] K. A. Meissner, Phys. Lett. B392, 298 (1997); hep-th/9610131 

[41] L. N. Granda, MPLA, DOI 10.1142/S0217732312500186, arXiv:1108.6236vl 
[hep-th] 

[42] L. N. Granda and E. Loaiza, Int. J. Mod. Phys. D21, 1250002 (2012), 
arXiv: 111 1.2454 [hep-th] 

[43] S. Capozziello, G. Lambiase, Gen. Rel. Grav. 31, 1005 (1999); gr-qc/9901051 

[44] I. P. Neupane, Class. Quant. Grav. 23, 7493 (2006); hep-th/0602097 

[45] T. Kobayashi, M. Yamaguchi, J. Yokoyama, Prog. Theor. Phys. 126, (2011) 511; 
arXiv:1105.5723 [hep-th] 

[46] A. De Felice, S. Tsujikawa, Phys. Rev. D84, 124029 (2011); arXiv: 1008.4236 
[hep-th] 

[47] A. De Felice, T. Kobayashi, S. Tsujikawa, Phys. Lett. B706, 23 (2011); 
arXiv: 1108.4242 [gr-qc] 

[48] A. De Felice, S. Tsujikawa, JCAP 1202, 007 (2012); arXiv: 11 10.3878 [gr-qc] 

[49] S. Nojiri, S. D. Odintsov, S. Tsujikawa, Phys. Rev. D71, 063004 (2005); hep- 
th/0501025 

[50] S. Nojiri and S. D. Odintsov, Phys. Rev. D70, 103522 (2004) [arXiv:hep- 
th/0408170]. 

[51] S. Cotsakis, I. Klaoudatou, J.Geom.Phys. 55 (2005) 306; arXiv:gr-qc/0409022. 

[52] S., I. Klaoudatou, J.Geom.Phys. 57 (2007) 1303, arXiv:gr-qc/0604029. 

[53] M. P. Dabrowski and T. Denkiewicz, AIP Conf. Proc. 1241, 561 (2010); 
arXiv:0910.0023[gr-qc] 



30 



